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I,  INTRODUCTION 


The  surface  waves  of  Rayleigh  were  first  studied  In 
detail  in  seismology  where  they  represent  a principal 
part  of  an  elastic  wave  coming  from  distant  earth  tremors. 
During  the  last  IO-I5  years  they  have  acquired  considerable 
significance  In  ultrasonics  for  flaw  detection  and  for 
locating  surface  defects. 

These  surface  waves  are  "attached"  to  the  boundary 
between  two  media,  one  of  which  must  be  an  elastic  solid, 
and  are  confined  to  a region  near  the  surface,  giving  little 
movement  at  great  depth.  These  are  important  at  great 
distances  from  the  point  of  origin  for  the  following  reason. 
A pulse  traveling  symmetrically  outwards  from  a point  source 
In  three  dimensions  gives  a displacement  Inversely 


proportional  to  the  distance  r from  the  source.  But  If  the 

disturbance  Is  confined  to  a given  depth  h,  (cylindrical 

symmetry)  Its  energy  over  a band  at  a distance  r Is 

distributed  over  an  area  Znrh  and  therefore  the  amplitude  of 

-1/2 

the  disturbance  will  vary  as  r . At  large  distances  the 

surface  wave  may  predominate  over  other  types  of  waves,  and 
this  actually  occurs  In  the  case  of  earthquakes. 

An  Interesting  application  of  surface  waves  is  In  the 


field  of  underwater  acoustic  wave  launchers  and  receivers.. 


Electromagnetic  surface  wave  antennas  are  now  used  In  radio  Sectibn  |0| 

Buff  Section 

communications  and  a soijiewhat  analogous  arrangement  may  be 

TION  

possible  In  sonar  work.  Present  day  acoustic  antennas 
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depend  upon  (excepting  linear  arrays)  their  cross  sectional 
dimension  to  control  the  beam  pattern.  To  obtain  a narrow 
pattern  a large  cross  section  (in  terms  of  wavelength)  is 
necessary.  ",  However,  large  frontal  areas  are  hydrodynamlcally 
incompatible  with  high  5>eed  underwater  vehicles. 3 In  a 
surface  wave  antenna  the  length  of  the  radiator  rather  than 
the  cross  section  is  of  primary  Importance  in  the  production 
of  a large  directivity  index.  But  before  any  worthwhile 
radiator  can  be  designed,  information  concerning  the 
transmission  of  surface  wave  energy  between  media  must  be 
available. 

Ihls  paper  deals  with  this  question  of  surface  wave 
transmission  specializing  in  the  field  of  launching  underwater 
acoustic  energy.  Included  will  be  a discussion  of  the 
mathematics  of  surface  waves  associated  with  boundaries 
between  a solid  and  a vacuum,  and  a solid  and  a liquid.  A 
description  is  given  of  the  experimental  setup  and  the 
resulting  beam  patterns  obtained  from  surface  waves  launched 
from  metallic  rectangular  bars  will  be  shown. 


II.  WAVES  IN  ELASTIC  MATERIAL 


Before  any  understanding  can  be  achieved  as  to  the 
nature  of  surface  waves  a general  understanding  of  the 
mathematics  of  waves  in  elastic  bodies  is  a necessity. 

A classic  reference  is  the  Theory  of  Elasticity  by  S. 

, (1)  A brief  review  of  waves  in  elastic  materials 

Timoshenko. 

will  now  be  undertaken. 

Assume  an  isotropic  cube  of  elastic  solid  as  shown 
in  Figure  1 on  the  following  page.  Summing  the  forces  in 
any  one  direction,  for  instance  the  x direction,  we  obtain 


£ F . 


da  ^ 6t  ^ 

+ — 6 '6  6 - a66  +(''■  + 6 ; 6 6 

X dx  xyyz  xyz  Vxy  dy  y y x z 


dT 


T 6 6 + ( T + 6 6 - T 6 6 

xy  X z V xz  Bz  x v xz  x 


y xz  X y 


(1) 


+ X6  6 6 - 0 

xyz 


where 


6 ,6  ,6 
x'  y ' z 


Normal  stress  parallel  to  x axis 
Stress  on  j face  parallel  to  i axis 
Component  of  body  force  per  unit  volume  in 
X direction 
Sides  of  the  cube 


Superscripts  refer  to  references  listed  in  the  bibliography 

> 
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Dividing  by  yields 


da  dT  dT 

^ ^ + X = 0 

dx  dy  oz 


(2) 


and  similarly,  summing  forces  in  the  Y and  Z directions 
yields 


dy  dx 


dT 


= 0 


dT 


= 0 


(3) 


(4) 


These  equations  must  be  satisfied  throughout  the  body 
including  the  boundary  surfaces.  Expressing  the  stresses 
in  terms  of  strains  yields  the  following  set  of  equations: 
(see  Appendix  A) 


uEX  . ^^x  _ , „ . 

J T * ~r7  “ A.i  + 2Ge 

X (u+l)(l-2u)  u+1  X 


(5) 


uEX 


Ee 


“ (u+l)(l-2u)  u+1 


\X  + 2Ge 


(6) 


Ec 

^EX  Z _ \ A op « 

(u+l)(l-2u)  u+1  z 


(7) 


1 


where 


V = Poisson's  ratio  = transverse 
contraction/elongation 
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= Uni  „ elongations  in  x,  y and  z directions 

X'  j z 

•C  = Volume  expansion  = e + e,,  + e„ 

X y z 

E = Modulus  of  elasticity  in  tension 
G = Modulus  of  rigidity  = E/(i;+l)2 
X = Lame's  constant  = tjE/(1+w)(1-2w) 

For  the  case  of  uniform  hydrostatic  pressure 
o = rr  = a = -P  and  ^ = e +p  -*-e_  which  is  equal  to 

X y z X y z 


r ^ 

va 

x 

w r ^ 

!!k. 

R 

R 

j » 1 ^ ^ ^ 

^ R R 

R 

)n-  - 

R 

R 

, l-2„ 

= -3(  ) P 

E 


We  can  now  solve  equations  (2),  (3)|  and  (4)  for 
particle  motion  by  using  stress  in  terms  of  strain  and  then 
strain  in  terms  of  displacement  as  given  below. 


P =:  ^ 

y 


(9) 


Substituting  equation  (5)  into  (2)  yields 

X -^  + 2G  ^ + g(^  + ^ ')  + ^ ) + X = 0 (11) 
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which  can  be  written  as 


>x  V J 


2 > > 

Where  7 ^ 

>y  ^2 


G7^u  + X = 0 


(12) 


Finally,  by  combining  terms  and  equating  this  to  particle 
motion,  the  equation  becomes 


(X+G)  ^ + G7^u  + X = 0 ^ 


>x 


St 


(13) 


Similarly 

(X+G)  If  + + Y = 0 ^ (14) 

(x+G)  If  + Gt^w  + Z = 0 ^ (15) 

^ St 

where  o = density. 

Now  assuming  no  body  forces  ( 1 . e. . X=Y=Z=0)  and 

v*7v4>v,by  differentiating  equation  (13)  with  respect  to  x, 

equation  (14)  with  respect  to  y,  and  equation  (15)  with 

respect  to  z,  and  adding  the  resulting  equations  we  can  write 
(2 ) 

In  vector  form 

0 A = (x+G)  71  + G7  •7A  (16) 

where  ^ = lu  + Jv  + kw 
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This  can  be  written  In  a more  favorable  form  by  using 
the  vector  Identity 

7 •7A  = -7X'TXA+7'’»A 

which  yields 


OA  = (x  + 2G)  77  ‘A  - G7  X X A (1?) 

This  Is  In  a favorable  form  for  the  examination  of 
two  special  cases;  (1)  A Is  an  irrotatlonal  vector, 
(7x7=0)  and  (2)  7 Is  a solenoldal  vector,  (7  . A « 0^ 

In  the  first  case,  where  A Is  Irrotatlonal,  equation  (1?) 
reduces  to 

0 A = (X  2G)  7^A  (18) 

This  says  that  a Is  propagated  as  a wave  with  velocity 
rj+ZG 

Vl  = V p • This  wave  Is  termed  a longitudinal  (or 
Irrltatlonal ) wave  since  the  only  propagated  disturbance  Is 
directed  along  the  direction  of  propagation  as  shown  In 
Appendix  B. 

In  the  second  case,  where  A Is  solenoldal,  equation 
(17 ) reduces  to 


’■*  ^ 2- 
PA  = G7  A 


(19) 


which  says  that  this  displacement  component  is  propagated  as 

/ 

a wave  with  velocity  = vG/o  . This  displacement  is 
directed  at  right  angles  to  the  direction  of  propagation  as 
is  shown  in  appendix  B and  is  known  as  a transverse  or  shear 
wave . 

It  is  Important  to  note  that  for  an  ideal  fluid  the 
shear  modulus,  G,  is  zero  and  therefore  a pure  shear  wave 
cannot  exist  and  a longitudinal  wave  would  have  a velocity 
of  = A/o  » 

Surface  Waves  in  Elastic  Material 

We  have  seen  that  an  elastic  solid  can  support  two 
types  of  waves J a longitudinal  and  a transverse  or  shear 
wave.  What  we  are  primarily  Interested  in,  however,  is  a 
wave  propagating  in  an  elastic  material  but  attached  to  the 
boundary  of  that  medium. 

In  his  fundamental  paper  on  surface  waves,  Lord  haylelgh 
concerned  himself  with  the  question  of  whether,  in  an  elastic, 
isotropic  seral-lnf inlte  solid  body,  a wave  system  could  exist 
which  decreased  exponentially  with  distance  to  the  boundary 
plane.  The  conclusion  he  arrived  at  was  that  such  a system 
is  possible  in  any  isotropic  medium  and  that  it  travels  in  a 
direction  parallel  to  the  boundary  with  a velocity  which  is 
smaller  than  the  velocity  of  transverse  waves. 

This  classic  paper  considered  a boundary  between  an 
elastic  solid  and  vacuum,  the  vacuum  being  Incapable  of 
supporting  any  acoustic  wave.  Assume  a longitudinal  wave 
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with  a velocity  potential  41  and  a transverse  wave  with  a 
velocity  potential  t,  in  the  elastic  solid  traveling  in  the 
X direction.  Using  the  plane  z=0  as  the  boundary  plane  with 
z taken  positive  toward  the  elastic  solid,  the  following 
boundary  conditions  must  be  satisfied  at  z=0: 

1)  Shear  stress  must  vanish,  thus 

If  * H)-  0 <20) 

2)  Normal  stress  must  vanish,  thus 

2Q  ^ = 0 (22) 


Writing  the  displacements  in  terras  of  4 and  y we 


obtain 


u = ^ + 31 


w 


>4) 

^2  " ’ 


0 


(23) 


Assuming  that  41  and  y are  in  the  form  of 

41  = A(z)  exp  [ik(x-ct)]  (24) 

'»  = B(z)  exp  [ik(x-ct)]  (25) 

where  k is  the  wave  number  2tt/>  and  substituting  equation 
(23)  into  equations  (13),  (14),  and  (15)  (the  equation  of 
motion)  the  following  conditions  result: 


A(z)  = C exp[-rz]  and  B(z)  = D exp[-sz]  (26) 


where  C and  D are  constants  and  r and  s are  the  positive 


11 


values  satisfying 


(27) 


where  a = 0 = V^,  and  c is  the  velocity  of  the  surface 

wave.  This  derivation  Is  shown  In  Appendix  C. 

Now  r and  s must  be  positive  real  numbers  In  order  to 
fulfill  our  definition  of  a surface  wave.  If  either  Is 
negative  or  imaginary  the  wave  amplitude  would  tend  to 
Infinity  with  depth.  A condition,  therefore,  for  the 

existence  of  a surface  wave  as  seen  from  equation  (2?)  Is 

2 2 2 
that  c shall  be  less  than  0 and  less  than  a . If  these 

conditions  are  not  satisfied  ♦ or  ? or  both  will  not 

represent  a surface  wave  but  a wave  traveling  upwards  or 

^ ^ (^) 
downwards . 

Using  equations  (23),  (24),  and  (25)  the  boundary 
condition  given  by  equation  (20)  now  becomes 


Since  In  an  elastic  solid  the  shear  modulus  0 is  not  equal 
to  zero,  we  can  write  this  boundary  condition  with  the  use 
of  equations  (24),  (25),  and  (26)  as 


-21krC|g|~^.^j  + D(s^  + k^)  = 0 


(29) 
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Equations  (30)  and  (32)  can  tlfien  de  combined,  as  Is  done  in 
Appendix  D,  to  yield  the  Eayleif^h  wave  equation 

2 r 6 4 2 o2  „2 

S r ^ 8 \ %(24  - 16  ^)  - 16(1  - %)1  = 0 (33) 

2 2 

This  is  a quadratic  in  terms  of  c /3  . One  of  the 
roots  is  obviously  equal  to  zero.  Hence,  from  equations  (2?) 
and  (30)  we  have  r=k=s,  C=-iD  and  then  substituting:  this 
into  equations  (24),  (25),  and  (26)  and  usins:  the  results  in 
equation  (23)  we  find  that  u=v=w=0.  Thus  there  is  no  wave 
at  all  and  this  is  not  acceptable. 


I 


( 1-  2» 
2(l-w 


Discarding  this  factor  and  writing 
the  Rayleigh  equation  now  becomes 
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8r3  - 


(»)] 


3^ 


- <1  - Irfer)  = “ 


(34) 


In  general  this  equation  has  one  real  root  and  two  complex 

roots  which  must  be  discarded  for  reasons  previously  discussed. 

Each  medium  has  Its  own  Poisson  ratio  which,  when  Inserted 

Into  equation  (3^),  will  determine  the  acceptable  value  for 
2 2 

c /3  . Once  this  ratio  is  determined  the  process  to  determine 
particle  motion  Is  straight  forward  and  Is  as  follows; 

1)  Determine  3^  from  3^  = G/o 

2)  Find  c from  the  solution  for  c^/3^ 

3)  Using  equation  (2?)  find  r and  s In  terms  of  k 

4)  Using  equation  (23)  find  u and  w In  terms  of  C,  D, 
k and  z 

Proceeding  In  this  manner  and  taking  the  real  part  of 
the  resulting  expression,  the  displacements  u and  w may  be 
written  as 

u = K(exp[-rz]  - M exp[-sz])  sin  k (x  - ct ) (35) 

w = K(N  exp[-rz]  - 0 exp[-ez])  cos  k (x  - ct)  (36) 

where  K Is  a constant  and  M,  N,  and  0 depend  upon  3,  C,  D, 

SCSI" 
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r and  So  These  may  also  be  written  as 


(5) 


0 ? 1/2  P P 1/2  1/2 

u = k(  11  .z-a^/aiL-  extC-kd  - oVb^)]  z 

^ (2  - c^/P^) 


1/2, 


- exp[-k(l  - c^/0^)]  ^j-sln(wt-kx)  (37 


w = -k| 


2(1  - c^/3^) 
2 - c^/3^ 


1/2 


? ? 1/2 

exp[-k(l  - c‘^/0‘^)  ]z 


1/2 


- (1  - c^/3^)  exp[-k (1-c^/P^  J z|  cos(wt  - kx)  (38 


In  Figure  2 on  the  following  page  are  plotted  the 
particle  orbits  In  a surface  wave  as  well  as  the  Poisson's 
ratio,  a/p,  and  c/3  with  which  they  are  associated.  Also 
with  reference  to  equation  (36),  the  depth  at  which  u=0 
Is  also  noted.  All  of  the  magnitudes  of  displacement  are 
referred  to  the  value  of  the  transverse  component  w at  the 
surface  of  the  solid.  As  Is  shown  In  Figure  2,  the  particles 
near  the  surface  travel  in  elliptical  orbits  In  a clockwise 
direction  as  viewed  In  a wave  moving  left  to  right;  at  a 
certain  depth  below  the  surface  the  particles  oscillate  In 
a vertical  line;  and  at  still  greater  depths  the  particles 
travel  In  elliptical  orbits  of  diminishing  dimensions  In 
which  the  direction  of  rotation  Is  counterclockwise. 

From  an  energy  standpoint,  about  95^  of  the  energy  Is 
contained  within  1 wavelength  of  the  surface. The 
experimental  results  of  Cook  and  VeinValkenberg  ^ are  shown 
In  Figure  3« 
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FIGURE  2 PARTICLE  ORBITS  IN  A 
SURFACE  WAVE 


Reprinted  from  Figure  3 of  "Surface  Waves  ot  Ultrasonic 
Frequencies"  by  S.  Cook  and  H,  Van  Van  Valkenberg  !-;lth  roT® 
additions . 
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+ + THEORETICAL  ASSUMES  V=-j 

00  EXPERIMENTAL  ALUMINUM 
25  MC 


DISTANCE  BELOW  SURFACE,  WAVELENGTHS 


FIGURE  3 

ENERGY  OP  SUBPAGE  WAVES  (PEB  CINT)  versus 
DISTANCE  BELOV  SURFACE  (WAVE  LENGTHS) 

Reprinted  from  Plgure  6 of  "Surface  Waves  at  Ultrasonlo 
Prequenoies”  by  1,  Cook  and  I.  VAn  Valkenberg 


III.  SUhFACE  WAVES  ON  A FLUID  LOADED  BOUNDARY 


In  the  preceedlnp  section  we  have  determined  the 
surface  wave  motion  for  a wave  ^tached  to  a boundary  between 
a solid  and  vacuum.  What  we  desire,  however,  Is  the  particle 
motion  for  a surface  wave  attached  to  the  boundary  between 
a solid  and  a liquid. 

The  procedure  here  Is  somewhat  different  than  that 
of  the  preceedlng  section.  This  is  broup'ht  about  as  an 
easy  way  to  handle  the  stress  at  the  boundary  which  are  now 
altered  by  the  appearance  of  a wave  in  the  liquid  medium. 

One  can  ®parate  the  wave  in  the  solid  Into  Its  longitudinal 
and  transverse  components  and  after  determlnlnp-  the  stress 
on  the  boundary  caused  by  each,  equate  them  to  produce  the 
proper  boundary  conditions.  This  will  then  yield  an  equation 
for  the  wave  traveling  along  the  surface  of  the  medium. 

In  an  ideal  fluid  the  modulus  of  rigidity,  G,  Is  equal 
to  zero.  Therefore,  the  fluid  can  support  a longitudinal 
wave  but  not  a shear  wave.  This  means  that  the  stress 
normal  to  the  boundary  will  be  continuous  while  the  stress 
parallel  to  the  boundary  must  vanish.  Using  the  coordinate 
system  as  shown  In  Figure  4 a longitudinal  wave  of  the  form 

X sln6.  + z COS0-1 

4i  = Aj^  exp[lv(t  - =)]  exerts  a stress  on 

the  plane  z=0  equal  to 

= -IvottA^  2 sln^e^  cotBj^  (39) 

= -IvottAj^  cos  20^ 


(40) 
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where  the  phase  velocities  have  been  omitted.  These  are 
obtained  from  equations  (20)  and  (22)  and  the  use  of 
Snell's  Law  sin  = 0/a  (sin  6^).  (41) 

In  the  liquid  medium  a wave 

X sin  6,'  + z cos  6,* 

((i'  = expLlv  (t  - ^)]  (42) 

appears,  which  Is  connected  with  the  movement  of  the  solid 
body  by  the  condition  that  the  vertical  motion  only  must 
be  continuous.  From  this  condition  It  follows  that 

slne^  sln6, 

= and  A,'  cos  6,'  = A,  cos  . (43) 

a*  a 1 i 1 1 

This  wave  In  the  fluid  exerts  on  the  Interface  a vertically 
directed  stress 

cotB-, 

'^z  = -ivo'a'A^  = -ivo’aA  (44) 

This  results  In  stress  differences  at  the  boundary  equal 
to 


*^zx  = 

-IvoCiA^  2 sln^G^  cos  6^ 

(45) 

o'  cot  6, 

0 

q 

N 

It 

1 

IvoaA,  (cos2e,'  - ) 

0 cot  ej^ 

(46) 

where  A now  denotes  Infinitesimal  change.  These  stress 
differences  can  be  represented  by  a vector  In  the  direction 
given  by 
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The  transverse  wave  in  the  solid  may  be  handled  in 
the  same  way.  Let  t be  of  the  form 
X sin  + z cos  6 

A2exp[lv(t  - ^ -)]  traveling  in  the  solid 

medium  and  incident  upon  the  boundary  as  shown  in  Figure  5» 
The  stresses  on  the  boundary  plane  z=0  caused  by  this  wave 
in  the  solid  are  by  use  of  equations  (20)  and  (22) 


(48) 

(49) 

The  accompanying  wave  in  the  fluid  is  propagated  in  the 
direction  6^  as  given  by 


= -ivo0A-  cos  26 
zx  2 2 

= +iVo6A2  2sin^  6^  cotd^ 


sin6^  a' 
sinS^  3 


(50) 


with  A^  satisfying  Aj^cose^  = -A2  sine2»  (51) 

This  wave  exerts  a stress  normal  to  the  boundary  equal  to 


ivo '3A2 

cote^ 


(52) 


The  stress  differences  then  caused  by  this  transverse  wave 
are 

ATzx  = -ivn3A2  cos  262  (53) 

“ Arr  = +lvo3A_(2sln^e  cot6  - — 2 ) (54) 

^ ocote^’ 
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The  stresses  due  to  the  transverse  wave  can  similarly  be 
represented  by  a vector  In  the  direction  given  by 


At 


zx 


cos  2t), 


Act, 


2sin  cotB^ 


oCOt©^ 


(55) 


As  the  normal  stress  on  the  Interface  mi’.st  be  continuous 
the  stress  differences  have  to  disappear  which  is  only 
possible  if  the  two  stress  differences  as  given  by  equations 
(47)  and  (55)  in  the  same  direction.  This  leads  to 


2sln‘^62COt6^ 


cos  2©, 


n 'cot©- 


( 56 ) 


cos  2© 


2 0 'cot©| 


2sin  e^oote^  - 


which  If  expanded  yields 

o h n'cot©-, 

cos'^  2©p  + 4sln^e„  cot©-,  cot©„  - = 0 (5?) 

^ z i z ocot©^ 


Although  this  condition  is  necessary  it  is  not 

sufficient.  The  stress  on  the  Interface  disappears  if  A^ 

and  A„  satisfy  the  equation  at  = a or 

^^long  ^^trans 

aA^2sln^y2COtaj^  = 0A2COs2©2  (58) 


Equation  (57)  determines  the  directions  of  the  three-wave 
system  possible  In  two  media,  one  of  which  Is  a liquid.  Such 
a system  will  be  a kind  of  Haylelgh  wave  If  the  waves 
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decrease  exponentially  with  increasing  distance  from  the 
interface.  In  this  case  cotB^^  and  cotB2  are  positive 
imaginary  and  cot 6^  is  negative  imaginary.  Now  by  use  of 
Snell's  law  and  knowing  that  c is  less  than  B,  equation  (5?) 


becomes 


c^) (a^ 


c2)l 


1/2^  ^ o'TlsLx-a?!  f 0 


This  is  the  dayleigh  wave  equation  as  expressed  in  equation 

(33)  but  in  a different  form  and  also  Includes  the  case 

where  the  second  medium  has  a non-zero  density. 

2 

One  root  is  again  c =0  which  as  before  is  not 

acceptable.  For  small  values  of  c^  the  left  hand  side  of 

equation  (59)  is  negative  and  if  c'"  is  taken  equal  zo  the 

smallest  of  the  three  velocities  (which  is  either  3 or  a' ) 

the  left  hand  side  will  now  be  positive.  There  will  therefore 

2 2 

always  be  a root  such  that  c will  be  less  than  a , 

For  small  values  of  n'a.'/na  the  transfer  of  eners-y 
is  also  small  and  propagation  will  not  be  much  different 
from  that  in  the  preceedlng  discussion  where  a vacuum 
replaced  the  liquid.  However,  in  the  case  where  a'  is  much 
smaller  than  3 and  also  smaller  than  c the  third  term  in 
equation  (59)  is  imaginary.  This  produces  a pair  of  complex 
roots  with  real  parts  equal  to  c. 

The  corresponding  wave  system  is  now  as  before,  but 
the  imaginary  part  of  the  roots  gives  rise  to  a slight 

CQIHMBM  DNIVERSITf 
K'!DS:N  ISBOSAIOillES 
CCiilACr  K'iinr-2BB(84) 
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exponential  decrease  with  distance  of  propagation  and  a small 

sinusoidal  variation  In  the  vertical  direction. ^ ^ ^ 

Another  factor  of  great  Importance  Is  that  a surface 

wave  can  only  be  propagated  In  a substance  whose  thickness 

( 7 ) 

Is  larger  than  two  surface  wavelengths.  Throughout  the 

preceedlng  sections  It  was  assumed  that  the  elastic  medium 
was  an  Infinite  half space,  but  In  any  experiment  or 
practical  case  this  cannot  occur.  As  the  thickness  of  the 
bar  Increases  the  distance  of  possible  surface  wave 
propagation  also  Increases.  For  a thickness  h = 3>j,  the 
distance  of  propap'atlon  Is  about  20>^,  for  h = 5>^  it  Is 
already  2000>^ 


f 


IV.  EXPERIMENTAL  ARRANGEMENT 


The  experiment  was  conducted  by  launching  surface 
waves  on  an  aluminum  bar  underwater  and  then  measuring 
the  pressure  pattern  in  the  liquid  at  the  end  of  the  bar. 

The  experimental  arrangement  shown  in  block  diagram  form 
in  Figure  6 consisted  of  two  main  parts;  the  acoustic 
equipment  and  the  electronics. 

The  steel  tank  in  which  the  exp'^rlment  was  carried 
out  is  6 1/4  feet  long  and  2 1/4  feet  in  width  and  depth. 

The  tank  is  lined  on  f'e  bottom,  ends  and  sides  with  two 
inch  thick  rubberized  hair.  Measurements  indicated  a 
normal  reflection  loss  of  26  db  for  the  material  at  a 

/ O \ 

frequency  of  1 mops'  . Since  a pulse  technique  was  used 
the  absorption  was  quite  sufficient  to  handle  all  possible 
reflections.  The  tank  was  filled  with  fresh  tap  water. 

As  was  asserted  earlier,  in  order  to  transmit  a 
significant  percentage  of  the  energy  into  the  fluid,  the 
ratio  should  be  as  large  a=^  possible.  With 

reference  to  Table  1,  giving  some  important  acoustical 
properties  of  metals  and  alloys,  it  is  seen  that  exceot 
for  bismuth  and  magnesium,  aluminum  appears  to  be  the 
best  choice  for  a favorable  amount  of  energy  transmission 
into  the  water.  Another  factor  dictating  this  choice  was 
that  bars  of  bismuth  or  magnesium  would  be  quite  exoenslve, 

A photograph  of  the  bars  used  as  well  as  pictures  of  the 
tank  and  electronic  equipment  are  shown  in  Figures  7,8  and  9, 
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TABLE  I 

PROPERTIES  OF  METALS  AND  ALLOYS 


Material 

Specific 

Gravity 

V 

long, 
ft , /sec. 

®H20 

° Met.  ^1 

Aluminum 

.3^ 

2.7 

16,740 

.11 

Alnlco 

.32 

7.0 

16,072 

.04 

Antlm.ony 

.33 

6.6 

11,152 

.07 

Beryllium 

.01 

1.8 

27,552 

.10 

Bismuth 

.35 

9.7 

2,624 

.20 

Brass 

.33 

7.7 

11,152 

.06 

Bronze  Phosphor 

.35 

8.8 

12,136 

.05 

Cadmium 

.30 

8.6 

*8  200 

.0? 

Cobalt 

.30 

8.7 

15,416 

.04 

Copper 

.35 

8.9 

11,670 

.05 

Duraluminum 

.33 

2.8 

16,400 

.11 

German  Silver 

.37 

8.1 

12,464 

.05 

Gold 

.42 

19.3 

5,717 

.04 

Iridium 

.33 

22.4 

4,920 

.05 

Iron 

.28 

7.8 

16,600 

.04 

Lead 

.45 

11.3 

15,700 

.03 

Magnesium 

.33 

1.7 

15,100 

.19 

Monel 

.32 

8.8 

14,760 

.04 

Nickel 

.31 

8.8 

16,072 

.04 

Palladium 

.39 

12,0 

10,496 

.04 

Platinum 

.39 

21,4 

8,815 

.03 

Rhodium 

.34 

l^,4 

16,072 

.03 

Silver 

.38 

10.5 

8,553 

,06 

Steel  C ,08 

.37 

7.7 

16,400 

.04 

Tantalum 

.31 

16.6 

11,152 

.03 

Tin 

.36 

7.2 

8,200 

.08 

Tungsten 

.17 

19.0 

14,104 

.02 

Zinc 

.43 

7.1 

12,146 

.06 
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Three  of  the  four  bars  used  were  machined  on  all  six 
faces  leaving  the  remaining  bar  to  compare  the  effect  of 
surface  roughness  on  surface  wave  propagation. 

The  transducer  used  to  launch  the  surface  waves  was 
a barium  tltanate  plate  In  a plexiglass  holder  which  was 
nounted  on  a plexiglass  wedge  as  shown  In  Figure  10.  The 
basic  equation  for  surface  wave  creation  Is,  with  reference 
to  Figure  10, 


, Q _ -1  ^longitudinal  wave  In  wedge  (60) 

sxn  t/,  — sin 

^surface  wave  In  bar 

and  therefore  V,  . . , . , , must  be  less 

longitudinal  wave  In  the  wedge 


than  V 


A sector  of  a circle 


surface  wave  in  the  solid* 
rather  than  a triangular  wedg^  with  a fixed  angle  was 
used  In  order  to  facilitate  the  use  of  different  metallic 
substances  if  this  was  desired.  Acoustical  contact 
between  the  transducer  housing,  wedge  and  bar  was  maintained 
by  a film  of  castor  oil.  The  loss  in  acoustic  pressure 
between  the  transducer  and  bar  (a  distance  of  4 inches) 
at  a frequency  of  760  kc,  which  was  used  In  the  experiment, 
is  approximately  20  db^^'^^  For  the  case  of  the  plexiglass 
wedge  having  a longitudinal  wave  velocity  of  9068  ft, /sec, 
and  aluminum  having  a surface  wave  velocity  of  9769 
ft. /sec,  the  optimum  Incidence  angle  Is  68°  which  was 
kept  constant  throughout  the  Investigation,  This  angle 
Is  the  solution  to  equation  (60)  and  Is  the  critical  wedge 
angle  at  which  maximum  conversion  of  a longitudinal  wave 


FIGURE  : R ' USE:.  .N  '"HE  EXPSRI>!EN1*' 


PIGUHE  8 EXPERIMENTAL  ARRANGEMENT 


FIGURE  9 ELECTRON I 


ARRANGEMENT 
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Into  a surface  wave  occurs^^^^.  This  transformation 

occurs  because  the  wedge,  acting  as  a prism,  sets  up  at 

the  boundary  a periodic  disturbance  with  a spatial  period 

of  > Since  the  incident  angle  is  greater  than  the 

angle  of  total  Internal  reflection  for  both  longitudinal 

and  shear  waves,  there  results  an  Inhomogeneous  wave 

known  as  a surface  wave.  The  wedge  method  was  also  used 

because  in  all  methods  of  surface  wave  excltetion  except 

the  wedge  method  the  resulting  wave  pulse  has  a longer 

( 12) 

duration  than  the  exciting  voltage  , 

The  receiving  probe  is  a Glennlte  Ultrasonic  Probe 
type  UP-SOO  manufactured  by  Gulton  Industries,  Inc,  of 
Metuchen,  N,  J,,  The  sensitive  element  of  this  unit  is 
a barium  tltanate  ceramic  cylinder  having  a diameter  and 
length  of  1/16  th  of  an  inch.  At  a frequency  of  76O  kc 
the  sensitivity  of  the  probe  is  -135  db  referred  to 
1 volt/dyne/square  centimeter. 

The  electronic  arrangement  is  Included  in  Figure  6 in 
block  diagram  form.  The  output  waveform  is  a pulsed  sine 
wave  with  a variable  carrier  frequency.  The  frequency 
and  amplitude  of  the  excitation  are  controlled  by  the 
Hewlett  Packard  Signal  Generator  model  606-A  [1]  while 
the  pulse  repetition  ratf  IS  determlrsd  by  the 
General  Radio  Unit  Pulser  model  1217-ASl  [2],  which  also 
acts  as  the  synchronizing  device  for  the  experiment.  The 


* Numbers  In  brackets  refer  to  circled  numbers  in  Figure  6, 
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frequency  rage  Is  50  kc  to  65  me  and  the  repetition  rate 
Is  adjustable  from  25  cps.  to  100k  cps . . The  pulse 
duration  Is  determined  by  the  General  Radio  Pulse,  Sweep 
and  Time  Delay  Generator  model  1391-B  [3],  with  £ variable 
pulse  width  from  0.025  usee,  to  100,000  usee,.  A gate  [4] 
whose  schematic  diagram  Is  shown  In  Figure  11  Is  used  to 
modulate  the  sine  wave.  After  gating,  the  pulsed  sine 
wave  Is  amplified  by  the  General  Radio  Power  Amplifier 
model  1233-A  [5]  with  a gain  of  60+  1/2  db  from  20  cps, 
to  2 meps,  and  Is  then  transmitted  by  the  barium  titanate 
transducer  [6],  The  output  wave  shape  has  a signal  to 
noise  ratio  of  500. 

The  signal  produced  by  the  receiving  probe  [?]  Is 
amplified  by  a specially  built  small  signal,  high  Input 
Impedance,  wideband  amplifier  [8]  with  a voltage  gain  of 
100  db.  The  schematic  of  the  amplifier,  which  has  a minimum 
signal  requirement  of  2 uvolt,  Is  shown  in  Figure  12, 

This  signal  Is  then  gatted  [9]  so  as  to  pick  out  only  that 
part  of  the  received  wave  produced  by  the  surface  wave. 
This  gate  Is  Identical  with  the  one  used  In  forming  the 
transmitted  signal,  but  Is  now  modulated  by  a square  wave 
produced  by  the  Rutherford  Electronics  Square  Wave 
Generator  model  B-2A  [lO],  This  generator  Is  triggered 
by  the  Pulse,  Sweep  and  Time  Delay  Generator  and  can  le 
adjusted  to  produce  a square  pulse  of  variable  width  and 
variable  time  delay  after  being  triggered.  The  resulting 
signal  Is  then  just  the  wave  produced  by  the  surface  wave 
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FIGURE  12  SCHEMATIC  OF  RECEIVING  AMPLIFIER 
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whan  the  variable  gate  is  adjusted  oorreotl/*  This  is 
then  detected  by  a peak  detector  £11] 'and'  then  recorded 
on  the  Oeneral  Badlo  Graphic  Level  Becorder  aodel  'S^Zlw^k 
Cl2]«  An  oscilloseope,  which  is  triggered  by  the  Uilt 
Pulser^  Is  used  to  view  the  waveshape  at  any  point  In  the 
systen.  This  Is  used  to  measure  pulse  width  and  amplitude, 
to  help  determine  wave  velocities  by  noting  time 
dt  'ferences  and  to  aid  in  positioning  the  gate  C9]«  Aie 
overall  sensitivity  of  the  receiving  system  was  such  as 
to  produce  an  output  of  ,1  volt  do/  80  mlorobar. 


V.  PROCEDURE 


In  order  to  obtain  the  greatest  sensitivity  it  was 
necessary  to  choose  an  optimum  frequency  and  select  the 
proper  Incidence  angle  for  the  transmitting  transducer. 

The  frequency  was  chosen  with  three  frequency  responses 
In  mind:  the  transmitter,  the  receiving  probe  and  the 
receiving  amplifier.  Also  taken  Into  consideration  was  the 
fact  that  the  smallest  dimension  of  the  bar  should  be  at 
least  5>  j,  to  Insure  propagation.  Since  the  aluminum  bar 
used  for  field  pattern  measurements  has  a surface  wave 
velocity  of  about  9800  ft, /sec,  and  a smallest  dimension  of 
1 1/4  Inches  the  frequency  chosen  must  be  greater  than  470  kc. 
The  frequency  chosen  was  760  ko.  At  this  frequency  the 
transmitter  had  a good  impedance  ratio,  the  probe  was  at  its 
greatest  sensitivity  and  the  amplifier  had  ample  gain. 

The  proper  angle  of  incidence  for  the  transmitter  is 
68°  as  given  by  equation  (60),  This  was  experimentally 
corroberated  and  the  wedge  was  positioned  18  Inches  from  the 
end  of  the  bar. 

Once  the  proper  settings  were  achieved,  care  was 
exercised  in  maintaining  them  throughout  the  experiment. 

Before  every  test  the  pulse  width,  pulse  amplitude,  gate 
width  and  power  supplies  were  checked.  In  the  experiment  a 
pulse  width  of  20  usee,,  a pulse  amplitude  of  35  volts  peak- 
to-peak  and  a repetition  rate  of  500  pulse  per  second  were 
used.  The  receiving  gate  width  was  set  at  25  usee,  to  insure 
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that  the  entire  wave  pulse  was  recorded  even  if  the  gate  was 
not  positioned  perfectly.  This  did  not  cause  any  trouble 
since  at  no  time  throughout  the  experiment  did  any  other  wave 
packet  arrive  at  the  probe  at  a time  less  than  5 usee,  from 
the  acoustic  wave  caused  by  the  surface  wave  In  the  bar. 

The  transmitter,  although  designed  to  produce  surface 
waves,  also  produced  two  other  observable  waves.  One  was  a 
longitudinal  wave  In  the  bar  and  the  second  was  a wave  In  the 
water  due  to  the  transmitter  configuration.  A shear  wave  In 
the  bar  was  also  produced,  but  since  all  otservatlons  were 
made  In  the  liquid,  this  was  not  apparent.  The  problem  now 
was  to  pick  out  that  part  of  the  wave  In  the  water  at  the 
end  of  the  bar  produced  by  the  surface  wave.  The  technique 
used  was  that  of  timing  the  various  wave  packets.  The  time 
of  arrival  of  each  wave  at  the  end  of  the  bar  could  be 
predicted  from  the  knowledge  of  the  different  wave  velocities. 
Once  the  waves  leave  the  end  of  the  bar  and  go  Into  the 
water  they  all  travel  with  the  same  velocity  and  preserve 
their  time  separation.  Since  the  longitudinal  wave  velocity 
In  aluminum  is  16,740  ft, /sec.  and  the  surface  wave  velocity 
In  aluminum  Is  9,800  ft. /sec.,  after  a propagation  distance 
of  18  Inches  the  leading  edges  of  these  two  wave  pulses  are 
60  usee,  apart.  The  water  wave  produced  by  the  transducer, 
traveling  at  5000  ft. /sec..  Is  now  210  usee,  behind  the  body 
wave  and  I50  usee,  behind  the  surface  wave  and  was  easily 
Identified.  Since  the  time  separation  betwen  wave  pulse® 
at  the  end  of  the  bar  was  at  least  60  usee,,  a pulse  width 


of  20  usee,  and  a gate  width  of  25  usee,  eould  be  used. 

These  predletlons  were  shown  to  be  oorreet  by  using  the 
oscilloscope  as  a time  measuring  device  and  at  the  end  of 
the  bar  all  wave  packets  could  be  clearly  distinguished. 

The  beam  pattern  measurements  were  taken  by  positioning 
the  probe  at  a given  height  (relative  to  the  top  of  the  bar) 
and  at  a given  distance  from  the  end  of  the  bar.  The  probe 
was  then  motor-driven  along  a path  normal  to  the  direction 
of  wave  propagation  moving  slowly  enough  so  that  the  strength 
of  one  received  pulse  was  not  significantly  different  from 
the  next  to  insure  proper  detecting  and  recording.  In  this 
manner,  by  changing  the  height  and  position  of  the  probe 
and  adjusting  the  gate  to  each  new  position,  the  acoustic 
field  produced  In  the  water  by  the  surface  wave  In  the  bar 
was  mapped. 
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VI.  RESULTS 


Experiments  were  conducted  using  4 different  bars.  All 
the  bars  were  of  aluminum  rolled  stock  and  were  the  following 
sizes:  1"  x 8"  x 2?"  machined,  1 1/4"  x 2"  x 2?"  machined, 

1 1/4"  X 2"  X 30"  unfinished  and  1/2"  x 1"  x 2?"  machined. 

See  Figure  7 for  a photograph  of  the  bars.  The  smallest  bar 
was  used  to  investigate  the  question  of  medium  thickness 
versus  distance  of  surface  wave  propagation.  The  two  bars 
with  a 1 1/4"  X 2"  cross  section  were  used  to  plot  the 
beam  pattern  as  well  as  to  observe  the  effect  of  surface 
roughness  on  surface  wave  propagation.  The  largest  bar  was 
used  to  see  what  effect  a large  medium  depth  had  on  the  beam 
pattern. 

The  pattern  produced  with  the  wedge  16  Inches  from  the 
end  of  the  1 1/4"  x 2"  x 27"  machined  bar  is  shown  in  Figures 
13,  14  and  15.  Figures  13  and  14  show  the  pattern  when  the 
1 1/4  inch  dimension  was  used  as  the  depth  and  Figure  15 
shows  the  pattern  when  the  2 inch  dimension  was  used  as  the 
depth.  The  surface  wave  launching  technique  used  also 
produced  a surface  vsave  on  the  bottom  of  the  bar  of  the  same 
duration  but  of  smaller  amplitude.  Since  surface  waves  can 

( 1 -3  ) 

travel  around  corners  , the  bottom  surface  wave  would 
turn  the  corner,  travel  across  the  end  of  the  bar,  and  inter- 
fere with  the  surface  wave  on  top  of  the  bar  before  the 
entire  top  wave  packet  had  reached  the  end  of  the  bar.  For 
this  reason  the  2 inch  dimension  was  then  used  as  the  depth 
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FIGURE  15  VERTICAL  CROSS  SECTION  OF  THE  RADIATION 
PATTERN  FROM  A 2"  THICK  BAR 


to  avoid  this  Interference  resulting  In  the  pattern  shown  In 
Figure  15.  Due  to  the  length  of  the  tank,  measurements  at  a 
distance  greater  than  20  Inches  from  the  end  of  the  bar  were 
unobtainable . 

Field  plots  showed  that  a beam  with  a width  approximately 
equal  to  the  width  of  the  bar  Is  produced.  This  Is  expected 
since  at  a frequency  of  760  kc  the  smallest  dimension  of  the 
bar  corresponds  to  8x^  while  the  2 Inch  dimension  corresponds 
to  13>j,  producing  a long  straight  line  source.  The  forked 
nature  of  this  pattern  as  shown  In  Figure  14  most  likely 
results  from  the  surface  wave  energy  being  concentrated 
along  the  edges  of  the  bar.  The  vertical  cross  section 
shown  In  Figure  I3  with  3 db  points  only  6°  apart  Is  not  so 
easily  explained,  nor  Is  the  angle  at  which  the  wave  Is 
propagated.  Since  almost  all  the  energy  of  a surface  wave 
Is  contained  within  one  wavelength  of  the  surface  the 
resulting  sharp  pattern  Is  surprising.  The  relatively 
small  (5°)  angle  of  the  main  lobe  relative  to  the  bar  axis 
raises  another  question.  The  answer  probably  lies  In  the 
fact  that  with  regard  to  equation  (59)  the  third  term  on  the 
right  hand  side  causes  the  particle  motion  as  shown  In 
Figure  2 to  be  tilted  with  respect  to  the  boundary.  This 
tilt  would  then  produce  the  observed  angle  of  propagation. 

The  pencil  like  nature  of  the  beam  has  two  plausible 
causes.  The  first  Is  that  since  there  Is  also  a surface 
wave  on  the  bottom  of  the  bar  and  these  wav«s  can  pi  opagate 
around  corners,  they  may  Interfere  with  each  other  producing 


the  observed  pattern.  Secondly,  the  beam  width  may  depend 

upon  wedge  distance  from  the  end  of  the  bar.  Thus  It  may  be 

possible  to  use  the  length  of  the  radiator  to  match  the  bar 

to  the  fluid.  In  his  paper  "Surface  Waves  In  Acoustics" 

(14) 

Brekhovsklkh'  states  that  concerning  the  case  of  surface 
radiation  from  the  end  of  a cylinder  that  by  Increasing  the 
length  of  the  radiator  It  Is  possible  to  make  the  width  of 
the  principle  lobe  of  the  directivity  characteristic  as  small 
as  we  like. 

A comparison  of  the  two  bars  of  the  same  dimensions, 
one  being  machined  on  all  6 faces  and  the  other  being  used 
unfinished,  showed  no  difference  in  the  resulting  pattern. 
However,  the  degree  of  attenuation  of  the  surface  wave  with 
distance  along  the  bar,  although  small  In  both  cases,  was 
greater  In  the  unfinished  bar.  This  was  observed  by  probing 
along  the  surface  of  the  bar  as  well  as  noting  the  difference 
In  the  strengths  of  the  resulting  fields. 

The  experiments  with  the  smallest  bar  revealed  that  the 
thickness  of  the  elastic  solid  Is  a definite  factor  In  the 
distance  of  propagation.  Using  the  1/2"  dimension  as  the 
depth,  which  corresponds  to  It  was  found  that  In  a 

distance  of  18  Inches  the  surface  wave  had  decreased  to  such 
an  extent  that  only  a very  weak  radiated  pattern  could  be 
obtained.  In  the  18  Inches  of  prooagatlon  , Its  strenth  had 
decreased  to  1/lOth  of  Its  original  value  and  a scalloped 
pattern  was  produced  due  to  the  size  of  the  bar  compared  to 
Jlp.  The  largest  bar  was  also  tried,  but  It  was  found  that 


although  a similar  pattern  as  shown  In  Figures  13,  14  and  15 
existed,  It  was  masked  by  the  radiation  produced  by  the 
surface  wave  traveling  across  the  end  of  the  bar. 


VII,  CONCLUSIONS 


The  conclusions  of  this  Investigation  are; 

1)  A surface  wave  can  be  propagated  along  a 
boundary  between  an  elastic  so'' Id  and  a liquid, 

2)  A pencil  like  beam  pattern  can  be  oroduced  by 
this  wave  when  launched  from  the  end  of  a 
rectangular  bar, 

3)  The  surface  condition  of  the  bar  Is  a factor  In 
determining  the  decrease  of  surface  wave  energy 
with  distance  of  propagation. 

4)  The  distance  of  orooagatlon  of  surface  waves  Is 
a fuctlon  of  the  thick  less  of  the  elastic  solid. 

There  still  remain  many  questions  which  warrant  further 
study.  These  Include  the  changes  In  beam  pattern  due  to 
variations  in  radiator  length,  radiator  dimensions  relative 
to  surface  wave  wavelength,  and  acoustic  Impedance  (ov) 
ratios  of  the  solid  and  the  liquid.  Up  to  now  many  more 
questions  have  been  raised  than  have  been  answered.  Some 
further  aspects  of  this  work  that  should  be  looked  Into  are: 

1)  Experimental  Investigation  of  the  directivity 
versus  radiator  length  by  measuring  the  field 
with  the  w“dge  at  different  distances  from  the 
end  of  the  bar, 

2)  An  expllnatlon  of  how  the  sui face  wave  prooagates 
around  a corner.  How  Is  the  particle  motion 
modified  by  the  discontinuity  In  the  elastic 


medium?  The  answer  to  this  could  shed  some  light 
on  the  observed  radiation  directivity, 

3)  By  use  of  some  other  suitable  materials  for  the 
elastic  sclld,  determine  the  effect  of  the 
acoustic  Impedance  ratios  on  energy  transmission 
into  the  fluid  and  beam  angle  of  propagation. 
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APPENDIX  A 


Obtaining  Stress  In  Terms  Of  Strain 


Using  Hooke's  law  stresses  can  be  written  In  terms  of 
strains  as  follows; 


= - »(a  or  = 


'’x  “ ”'’y  - ’'"z 


(A-1) 


= l['’y  “ ^«y  = “ 


”‘"x  •’y  - 


(A-2) 


•z  “ e["z  " * a 


(A-3) 


Solving  for  & 


E«. 


Ee. 


- tf  - 1> 


1 -1 


= 


-V 

-u 


• w 
1 


- V 

- 

1 


Ee  + u^Ee  + ij^Ee  + vE* 

A y z z 

1 . „3  . „3  . „2  _ „2 


•’x  - 


-2,,^  - 3’’^  + 1 


= 


e|^»^(1-i>^)  + i>(l-»-v)i>^  ♦ 
(»-H)(w+l)(-2w-H) 


(A-i+) 
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now  = (l+u ) (l-2w ) + vil'^v)  , Therefore 


- 


Er('i+„)(i-2»)  + + ffd+v)  («„+«,)'' 

— .1.  i iL ^ ^ 

(»+l)(„+l)(l-2n) 


(w*-l) 


■*- 


(u+l)(l-2-0 


= ^ fix 

^ (u+l)(l-2v)  »+l 


tjE,.  Ee 

Similarly  o = + — ^ 

y (w+l)(l-2„)  ,,+1 


= £ + — £ 

^ (i,+1)(1-2t,)  w+1 


(A-5) 

{A-6) 

(A-7) 

(A-8) 
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Displacement  Directions  Of  Lon.^^ltudlnal  and  Transverse  Waves 

With  reference  to  equations  (18)  and  (19),  the 

direction  of  particle  displacement  can  be  found  for  both 

longitudinal  and  transverse  waves. 

I? ) 

a)  Longitudinal' 


which  means  that  u,  v,  w are  propagated  with  the  same  velocity 
( Irrotatlonal ) . Assume  u,  v,  w are  functions  of  x and  t only 
and  do  not  depend  upon  y and  z.  Equation  (B-1)  then  reduces 
to  the  three  equations 


..  V ^^v 


w j- 


(B-2) 


now  since  ^ is  irrotatlonal  '’x^=0  and 


vw  _ ^ ^ - iii 

^y  ■^z  az  ~ ^x  >x  ~ >y 


but  since  u,  v,  and  w are  functions  of  x and  t only 


lx.  = _ 2ii  = ^ 

>x  ■'x  ~ ^y  ^z 


0 and  hence  v Is  not  a function  of  x and 


u Is  a function  of  x oniy. 

Therefore,  the  only  propap-ated  displacement  Is  u,  which 
Is  directed  alona  the  direction  of  propagation. 
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( 2 ) 

b)  Transverse'  ‘ 


Assuming  as  above  that  a depends  on  x and  t alone  and 
since  now  a Is  solenoldal (7* a=0 ) we  readily  see  that  here 

= 0 and  hence  u Is  not  propagated  as  a waveo  This  means 

o X 

that  the  propagated  displacement  Is  perpendicular  to  the 
direction  of  propagation. 
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Evaluation  Of  A(z)  and  B(z) 

In  order  to  find  the  exact  form  of  the  velocity 
potentials  4)  and  y we  must  determine  A(z)  and  B(z). 

Substituting  Equation  (?3)  into  Equation  {I3)  with 
no  body  forces 


(c-1) 


now 


?>  ^ _ > ^^u  ^ ^ ^ aw~N  _ ^^u  ^ ^^v  a 

^x  ^y  >zy  ^^2  >x>y 


w 


^ ^ 0 * -iifii.  - 

^^2  V >X  ?Z  ■ ^X^Z  ^Z  >X 


Ji/ih  + 

^ :»x2  :.z2  ■ 


■'x'zy 


(C-2) 


• > Zj. 

^ 


(C-3) 


Separating  4)  in  v in  (C-1) 


dt 


2 


= (>+G)  7^41  + 

^x  ^ ' ^xy 


(C-4) 


n I-  = (•v+G)'’^4i  + G'’^4>  = (>+2G)'’^4) 

dt'^ 


(C-5) 


» ^ Cl?  = 


> '»  „ 2 
n — p = G" 


Substituting  (26  A)  Into  {C-6) 


(C-6) 


^ 2 

= -ok^c^A{z)  exp  ik{x-ct)  = - A(z  )k^exp  ik(x-ct ) 


- k^c^A(z)  = - k^A(zO 


(C-7) 


where  a = 7 = V, 

' 0 b 


A(z)  = C exp[-rz] 


where  G is  a constant. 

Wow  substituting  (26  B)  into  (C-6) 


(C-8) 


^ 2 
n = -nk^c^B(z)  exp  Ik(x-ct)  = G(^-k^B(z)  + Ik(x-ct) 


■ k^c^B(z)  = ay^g-Lz),  _ 1^23(2)^ 

dz"^ 


B ( z ) = D exp[-sz] 


(C-9) 


(C-10) 


2ikrC 

2*,, 2 
s +k 


where  e = = Jo/n 


and  D 


(C-11) 
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Substituting  (C<-8)  into  (C-7)  yields 

-k^c^C  expC*rs]  * a^(r^C  expC-rz3  • k^C  exp[»rz'H 

^ = 1 - c^/a^  <2? A) 

k* 

and  substituting  (G«10}  into  (C«9>  yields 

-k^o^D  exp[-sz]  = 8^{s^D  expC-sz]  - k^D  expC-s*3> 

2 

^ « 1 - 0^/0^  (27B) 


APPENDIX  D 


57 


Derivation  Of  The  Rayleigh  Wave  Equation 

Knowing  the  boundary  conditions,  the  velocities  of  the 
longitudinal  and  transverse  waves  and  the  displacements  in 
terms  of  the  velocity  potentials  the  Rayleigh  wave  equation 
can  be  derived. 

The  velocity  equations  can  be  written  as 


> + 20  = pa^,  G = 

Using  equation  (C-11)  which  is  D = ^ 

8 

and  substituting  It  Into  equation  (32)  yields 
+ (>+2G)r^C  + = 0 


or 

r2 

Ers 

■ G * ^ nr-J 

F; 

1 

From  (D-1) 

> = 

iia^-20 

•4  - S - 

2 = 

ai.  2 

and 

>-»-2G 

32 

(D-1) 

(D-2) 


(D-3) 


(D-4) 


(D-5) 
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Using:  the  compatibility  equations  (27)  and  substituting  (D-4) 

I . 

and  (D-5)  Into  (D-3)  yields 


2 


c_  _ 4rs 
0^  s^+k^ 


(D-6) 


Squaring  both  sides  and  using  equation  (2?)  for  and  s^ 
yields 


l6k^(l-cVa^)(l>o^/B^) 
k^  ( 4-40^3  ^ cVb^) 

^ l6(l-c^/3^){l-cV0^) 


16(l-oW)(l-c^/B^) 

(2-cVb^)^ 


(D^7) 


Expanding  both  sides  and  setting  the  right  hand  side  equal 
to  zero  results  In  the  Haylelgh  wav*  equation  which  Is 


c 


16  . i6ri . = 


0 


(33) 


